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Abstract 

Since the definition of the gluon to hadrons fragmentation amplitude < H,X\g > involves 
gluon field Q^°'{x) corresponding to single gluon incoming state \g >, the present dehnition of 
the gluon fragmentation function at high energy colliders [which involves non-abelian field tensor 
piJ-i'a^x) = d'^Q’^°‘{x) — d'^Q^°‘{x) + gf°‘^^Q^^{x)Q^^{x) instead of Q^°‘{x) in the initial state] is 
not consistent with the single gluon incoming state \g >. In this paper we derive the correct 
definition of the gluon fragmentation function at high energy colliders from first principles which 
is obtained from the single gluon incoming state \g > and is gauge invariant and is consistent with 
the factorization theorem in QCD. 
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I. INTRODUCTION 


In order to study hadron production at high energy colliders one needs to know how a 
hadron is formed from partons. However, this is an unsolved problem in QCD so far because 
the formation of hadron from partons involves non-perturbative QCD which is not solved 
yet. Hence, at present, we depend on experimental data to study hadronization in QCD. 

The usual procedure in the high energy phenomenology is to determine the parton to 
hadron fragmentation function from a set of experiments at some momentum transfer scale 
and then determine its evolution to other momentum transfer scale by using DGLAP evo¬ 
lution equations |^. Once the fragmentation function is determined in this way then one 
uses it to predict physical observable at other collider experiments. Hence in order for this 
prescription to work it is necessary that these fragmentation functions are universal, i. e., 
they do not change from one collider experiment to another collider experiment. From this 
point of view it is necessary to use the correct dehnition of the fragmentation function at 
high energy colliders. If one does not use the correct dehnition of the fragmentation function 
at high energy colliders then one will obtain divergent cross section for physical observable 
at high energy colliders. This can be seen as follows. 

Consider, for example, the hadron production at high energy colliders. If the factorization 
theorem holds 5| then the formula for the hadron production cross section at high energy 
hadronic colliders is given by 


a{AB^H + X) = Y. 

a,b 


dxi j dx2 J dz fa/Aixi^Q"^) fb/B{x2,Q‘^) d(a 6 -)■ cd) Dh/c{z,Q‘^ 


( 1 ) 


where fa/A{.Xi,Q‘^) (or fb/B{x2,Q‘^)) is the parton distribution function (PDF), dl^ab -A- cd) 
is the partonic level scattering cross section and Dh/c{z, Q"^) is the fragmentation function 
(FF) where a, b,c,d = q, q, g represent quark, antiquark, gluon. 

Hence if one does not use the correct dehnition of the fragmentation function in eq. 
o then any uncanceled infrared divergences in the Feynman diagrams in the partonic 
level scattering cross section d{ab cd) can not be absorbed in the dehnition of the non- 
perturbative fragmentation function DH/ai^iQ"^) in eq. ([I]). Hence it is necessary to derive 
the correct dehnition of the fragmentation function Dfj/a{z,Q‘^) from hrst principles. 

A correct dehnition of a parton to hadron fragmentation function at high energy colliders 
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must satisfy the following three properties: 1 ) since we are considering a single parton 
fragmenting to hadrons, it must be obtained from the single parton state |c >, 2 ) it must 
be gauge invariant and 3) it must be consistent with factorization theorem in QCD. The 
present dehnition of the quark to hadron fragmentation function at high energy colliders is 
consistent with these three properties However, the present definition of the gluon to 
hadron fragmentation function at high energy colliders is not consistent with these three 
properties as we will see below. 

The present dehnition of the gluon to hadron fragmentation function at high energy 
colliders is given by 

fOO 

[<0|F+^“(0) X PexpHr?T(^)' / dz-Q+\0,z-,0T)]aUP^,0T)aH{P^,0T) 

Jo 

roo 

X Vexp[-tgT^^^^ dz-Q+'=(0, y" + z", Ot)] x F^+“(0, i/", 0 t)|0 >] (2) 

where 

F‘Ax) = 9Mx) - a,Q‘{x) + gr'’Vi{x)Qi{x), (3) 

Pbc'^’^ = and is the quantum gluon held. 

The expression Vexp[igT^^'>^ dz~Q~^^{0, z~,0 t)] in eq. ([ 2 ]) stands for the light-like 

^1 

Wilson line. Since the light-like Wilson line in QCD involves soft-collinear gluon held 
it can be treated classically as y4^“(a;) becomes the SU(3) pure gauge when Wil¬ 
son line becomes light-like (see section III for details). Hence for the better under¬ 
standing of the single gluon initial state > in eq. ([ 2 ]) let us replace the light-like 

Wilson line Vexp[igT^^'>^ dz~Q~^’^{0, z~ ,0t)] in eq. ([ 2 ]) by its classical counterpart 

^1 

VexpligT^^^'^ dz~z~,0 t)] where is the classical SU(3) pure gauge back- 

ground held. With this one hnds from eq. ([3]) that in eq. ([2]) contains 

quadratic, cubic and quartic powers of the (quantum) gluon held Q^°‘{x) which implies that 
the present dehnition of gluon fragmentation function at high energy colliders in eq. ([ 2 ]) is 
not consistent with the single gluon initial state l^f >, he., F'^^“(x)F“j^( 0 ) 9 ^ in 

eq. ([2]) in QCD. This implies that the present dehnition of the gluon fragmentation function 
at high energy colliders is not correct. 

In this paper we derive correct dehnition of the gluon fragmentation function at high 
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energy colliders from first principles which is obtained from the single gluon initial state 
\g > and is gauge invariant and is consistent with the factorization theorem in QCD. 

We hnd that the correct dehnition of the gluon fragmentation function at high energy 
colliders obtained from the single gluon initial state \g > and is gauge invariant and is 
consistent with the factorization theorem in QCD is given by 

zk'^ 


D 


H/g 




poo 

[< 0|(5^“(a:“) X Pexp[—/ dz~+ z~)] ,0T)aH{P'^,0 t) 

Jo 

X Vexp[-igT^^'>^ r dz-A+%z-)] X Q“(0)|0 >] (4) 

which is valid in covariant gauge, in light-cone gauge, in general axial gauges, in general 
non-covariant gauges and in general Coulomb gauge etc. respectively. Note that the correct 
dehnition of the gluon fragmentation function at high energy colliders in eq. (jl]) is gauge 
invariant with respect to the gauge transformation 


TM“(x) = U(x)r‘Al(x)U-'(x) + -ldM^)]U-'{x), 

'^9 


u{x) = (5) 


along with the homogeneous transformation 


BJa 




Uix) = 


( 6 ) 


where Q^°‘{x) is the hard (quantum) gluon held which fragments to hadron PI and A^°'{x) 
is the soft-collinear gluon held which is the SU(3) pure gauge background held. After 
renormalization the gluon fragmentation function is expected to obey a QCD evolution 
equation, like DGLAP equation Ij, which follows from renormalization group equation. 

We will provide a derivation of eq. (jl]) in this paper. 

The paper is organized as follows. In section II we discuss the held theoretical derivation 
of the dehnition of the gauge non-invariant gluon fragmentation function in QCD. In section 
III we discuss the gauge invariant non-perturbative gluon correlation function and proof of 
factorization theorem in QCD. In section IV we derive the correct dehnition of the gluon 
fragmentation function at high energy colliders as given by eq. (|1|). Section V contains 
conclusions. 
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II. DEFINITION OF THE GAUGE NON-INVARIANT GLUON FRAGMENTA¬ 
TION FUNCTION IN QCD 

In this section we will present the derivation of the definition of the gange non-invariant 
gluon to hadron fragmentation function in QCD. The derivation of the dehnition of the gauge 
invariant gluon to hadron fragmentation function at high energy colliders which is consistent 
with factorization theorem in QCD will be discussed in next sections. For application to 
collider experiments it is convenient to use light cone quantization formalism. To introduce 
the concept of derivation of the dehnition of the fragmentation function using quantum held 
theory we will consider the scalar gluon hrst before considering the gluon in QCD. Note that 
the scalar gluon is not physical. The only reason we have considered it here is for simplicity 
in order to motivate the similar derivation of the dehnition of the gauge non-invariant gluon 
to hadron fragmentation function in QCD. 


A. Definition of Scalar Glnon Fragmentation Function Using Quantum Field The¬ 
ory 

Consider a scalar gluon with four momentum in vacuum fragmenting to a hadron 
with four momentum P^. Note that the derivation of the dehnition of the scalar gluon 
fragmentation function is described in js, 111 which we will briehy review here in order to 
motivate the similar derivation of the dehnition of the gauge non-invariant gluon to hadron 
fragmentation function in QCD. In the light-cone coordinate system the scalar gluon held 
(j){x~,XT) can be written as 


(j){x ,xt) = 


dk+ 


(27r)3 J yj2k+ 


d^kr kr) + kr)] 


( 7 ) 


where a is the annihilation operator and is the creation operator. The single particle 
scalar gluon initial state \k~^,kT > can be obtained from the vacuum state |0 > by using 


\k~^, kT >= a^kd^^ A;t)|0 > 


( 8 ) 


where 


a{k'^, A:r)|0 >= 0. 


(9) 
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The normalization condition is given by 


< fc+, k'T >= {2'Kf5{k'^ - k'^)5^^\kT - k'r). (10) 

As explained in [h| the correct interpretation of the state \k'^^kx > is created by an appro¬ 
priate Fourier transform of the corresponding held operator and should not be associated 
with on-shell condition /c^ = 0 of the gluon. 

For the inclusive production of a hadron H created in the owf—state \H,X > from the 
scalar gluon in the m—state \k'^,kT > the probability amplitude is given by 


< H,X\k+,kT > 


( 11 ) 


which gives the probability distribution hk{P) of the hadron with momentum P from the 
parton of momentum k 

hk{P) < kT\k'+, k!r>=Y. < ^t\H, X >< H, X|fc'+, k'^ > 


X 


= Y. < fc+, fcr|ayP)|X >< X\aH{P)\k'^, k'^ >=< fc+, fcr|ayP)aH(P)|fc'+, k'^ > 

X 

( 12 ) 

where is the creation operator of the hadron H. Hence we hnd (by using eqs. ([H]) and 
ffT0|) l that 

h{z, Pt) < A;+, krW^, k'r >= 2z{2xfDHia{z, PT){2xf6{k-^ - k'^)5^^\kT - k'r) 

=< ^\a{k^,kT)a}H{P^,PT)aH{P^,PT)a\k'^,k'T)\^ > ( 13 ) 

where Dn/aiz, Pt) is the transverse momentum dependent fragmentation function and = 
is the longitudinal momentum fraction of hadron H with respect to parton a. Using 


{271)''^ < 0\(j){x ,a;r)al^(P+,Pr)a//(P+,PT)0(O)|O 
< 0|e-*^^^'+*^^'^^a(fc+,fcT)ayP+,PT)aH(P+,PT)a^(fc'+,fcT)|0 > 


dk~^ 


■.(fki 


dk 


!+ 




d^k'l 


( 14 ) 


we hnd from eq. (B 

DhiAzXt) = U< 0|.^(x-,XT)apP+,Pr)«H(i’+.Pr)^(0)|0 > . (15) 

We rewrite the dehnition in a form analogous to the dehnition of the distribution of partons 
in a hadron where the transverse momentum is of the parton relative to the hadron rather 
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than vice versa. Hence making a Lorentz transformation to a frame where the hadron’s 
transverse momentnm is zero: 


(F+,Pt) = (P+,0) 

(fc+,0) = (fc+,-PT/^). (16) 


we hnd from eq. (na 

DH/a{z,PT) = —] < O|0(a:-,a:T)aL(P+,OT)aH(P+,OT)0(O)|O > . 

Integrating over pt we hnd that the scalar glnon to hadron fragmentation fnnction is given 
by 

Dn/aiz) = ^ /< O|0(a;-)al,(P+,OT)aH(P+,OT)0(O)|O > . (18) 

B. Definition of Gauge Non-Invariant Gluon Fragmentation Function in QCD 

In analogons to eq. ffTSD for the scalar glnon fragmentation fnnction we hnd that the 
dehnition of the gange non-invariant glnon to hadron fragmentation fnnction in QCD is 
given by 

DH/g{z) = ^ /da;-p^^"-/^<0|Q^“(a;-)aUP+,0T)ar,(P+,0r)Q“(0)|0> (19) 

which contains Q^“(a:)(5))(0) which is proportional to the qnadratic power of the (qnantnm) 
glnon held instead of the non-abelian held tensor P^^“(a:)P“j,(0) as in eq. ([2]) which 

contains qnadratic, cubic and quartic powers of the (quantum) gluon held where F^^{x) is 
given by eq. ([3]). Hence we hnd that, unlike eq. ([2]), the gauge non-invariant dehnition of 
the gluon to hadron fragmentation function in eq. flT^ is obtained from the single gluon 
initial state > in QCD. 

In order to obtain gauge invariant dehnition of the gluon to hadron fragmentation function 
which is consistent with factorization theorem in QCD we need to incorporate Wilson lines 
appropriately in the dehnition of gauge non-invariant gluon to hadron fragmentation function 
in QCD in eq. flT^ . Hence we need to prove factorization theorem of the non-perturbative 
gluon correlation function in QCD. 


( 17 ) 
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III. GAUGE INVARIANT NON-PERTURBATIVE GLUON CORRELATION 


FUNCTION AND PROOF OF FACTORIZATION THEOREM IN QCD 


As mentioned above the definition of the gluon fragmentation function in eq. fll9p is 
not gauge invariant. Hence it can not cancel any uncanceled infrared divergences in the 
perturbative cross section a{ab —)■ cd) in eq. ([I]). This implies that we need to prove 
factorization theorem for the infrared divergences of the non-perturbative gluon correlation 
function of the type < 0 |Q“(xi)Qy(x 2)|0 > in QCD in order to obtain the gauge invariant 
definition of the gluon fragmentation function from eq. 0191) . 

The infrared divergences can be studied by using eikonal approximation in quantum field 
theory. The Eikonal propagator times the Eikonal vertex for a gluon with four momentum 
interacting with a light-like quark moving with four velocity is given by 


12 


22 | 




( 20 ) 


I ■ k + ie 

Hence for a light-like quark attached to infinite number of gluons we find the eikonal factor 
d^k I ■ A^{k) 2 r d^k2 TH ■ A%ki)TH ■ A\k2) 


l + grj 


te J 


+ ... 


(27r)^ I ■ k + ie J (27r)"^ (2vr)"^ {I ■ ki + ie){l ■ {ki -|- ^ 2 ) + ie) 

roo POO POO 

= 1 + igT'^ / dXl ■ A\l\) + ^2^2 / / d\2T^l • A^{1\i)TH • A\l\2) + ... 

^0 JO JAi 

= 1 + igT’^ / d\l ■ A\l\) + / d\i / d\2Tn ■ A‘^{1 \i)TH ■ A\1\2) + ... 

7o 2! Jo Jo 

POO 

= V exph^ / dXL ■ A“(/A)T“] (21) 

Jo 


which exactly describes the infrared divergences arising from the infinite number of soft 
gluons exchange with the light-like quark where V is the path ordering and the gluon field 
A'^°‘{x) and its Fourier transform A^^°‘{k) are related by 

d^k 


A^\x) = 


( 2 ir) 




ik-x 


( 22 ) 


. A light-like quark traveling with light-like four-velocity d produces SU(3) pure gauge 
potential A^°‘{x) at all the time-space position x^ except at the position x perpendicular to 
the direction of motion of the quark (/ • T = 0 ) at the time of closest approach 12 , 24 1 . 

When A>^°‘{x) = as in eq. fl^ we find l-x = Xl-l = X^0 which implies that the 

light-like quark finds the gluon field in eq. (1^ as the SU(3) pure gauge. The SU(3) 

pure gauge is given by 


T“A“(x) = -[d,U{x)] U-\x), 
'^9 


U(x) = 


(23) 














which gives 


U{xj) = 


(24) 


Hence when A^°‘{x) = A^°‘{Xl) as in eq. fl2lD we find 


line in QCD for infrared divergences is given by 


igsl 


rom eq. 


^ p^-ig d\l-A-(x+lX)T- ^ pgT-u,-{x)_ 


that the light-like Wilson 


(25) 


In the adjoint representation of SU(3) the corresponding path ordered exponential is given 
by 

roo 

d\l-A\x + l\)T^^'^% (26) 

Jo 

which under non-abelian gauge transformation, as given by eq. ([5]), transforms as 

$l,(a;) = H dXl ■ A^{x + IX)T^^^% U^dj^x) = (27) 

^ 0 

Hence the effect of infrared gluons interaction between the partons and the light-like Wilson 
line in QCD can be studied by putting the partons in the SU(3) pure gauge background field 
which implies that the infrared behavior of the non-perturbative gluon correlation function 
of the type < 0\Q^{x)Q^^^{x')Qx{x")...\0 > in QCD due to the presence of light-like Wilson 
line in QCD can be studied by using the path integral method of the QCD in the presence 
of SU(3) pure gauge background field as given by eq. 


The non-perturbative gluon correlation function of the type < 0\Q‘l^{xi)Ql{x2)\0 > in 


QCD is given by 


,Q,C 


< 0|Q“(a;i)Q^(a;2)|0 >= 


QI{xM{x2) 


X 


^ ’ 


(28) 


where the suppression of the normalization factor Z[0] is understood as it will cancel in the 
final result [see eq. fl3^ ]. 

Background field method of QCD was originally formn 
extended by Klueberg-Stern and Zuber j^, and by Abbott 
which can be useful to construct gauge invariant non-perturbative green’s functions in QCD. 
This formalism is also useful to study quark and gluon production from classical chromo 


ated by ’t Hooft and later 
6j . This is an elegant formalism 
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field 


30| via Schwinger mechanism [^1, to compute (3 function in QCD 32|], to perform 


calculations in lattice gauge theories |33| and to study evolution of QCD coupling constant 
in the presence of chromoheld 341. 


The non-perturbative gluon correlation function o 


the background held method of QCD is given by 


je-g 


^e type < 0 \Ql{xi)Ql{x 2)\0 >a in 


13| 


< 0|Q“(a;i)Q^(a;2)|0 >a= J [dQ][#][#] Q‘I,{xi)QI{x2) 


X 


5oj’^ 


(29) 


where the suppression of the normalization factor Z[0] is understood as it will cancel in the 
hnal result [see eq. (13^ ]. In eq. (j^ the gauge hxing term is given by 

(30) 


G“(0) = = D^[A]QI“‘ 

which depends on the background held A^°‘{x) and 

KAA +01 = d,[Ai + Qi\ - s„|2i» + ca+gr>-[Ai + oy K + o:i 


(31) 


When the background held is the SU(3) pure gauge as given by eq. (123|1 then we hnd 


that Uh 


28| 


< 0\Q‘(xi)Qt(x2)\0 >=< 0\‘t^Axi)Ql(xi)^id(x2)Qi{x2)\0 >A 


(32) 


which proves factorization of infrared divergences at all order in coupling constant in QCD 
where is the non-abelian phase factor or eikonal factor in the adjoint representa¬ 

tion of SU(3) as given by eq. fl26l) . Note that eq. fl3^ is valid in covariant gauge, 
in light-cone gauge, in general axial gauges, in general non-covariant gauges and in gen¬ 


eral Coulomb gauge etc. respectively [26 


From eqs. 


and (j2711 we hnd that 


< 0 \^ac{.xi)Q‘l^{xi)^bd{x 2 )Qt{x 2)\0 >A in eq. fl3^ is gauge invariant and eq. fl3^ is con¬ 
sistent with the factorization of infrared divergences at all order in coupling constant in 
QCD. 


IV. CORRECT DEFINITION OF THE GLUON FRAGMENTATION FUNCTION 
AT HIGH ENERGY COLLIDERS 

From eqs. (fHj), (13^ and (1261) we hnd that the correct dehnition of the gluon fragmentation 
function at high energy colliders obtained from the single gluon initial state \g > and is gauge 
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invariant and is consistent with the factorization theorem in QCD is given by 



[< 0|Q'^“(x)PexpH^ r dXl ■ A\x + , 0T)aj^(P+, Ot) 


poo 

Vexp[-tg dXl ■ A^(/A)T(^)'=]Q“(0)|0 >].+=.,=o 


X'^=Xj'=0 


( 33 ) 


which is valid in covariant gauge, in light-cone gauge, in general axial gauges, in general 
non-covariant gauges and in general Coulomb gauge etc. respectively. Eq. (13^ can be 
written as 



(34) 


which reproduces eq. (jl]). 

This completes the derivation of the correct dehnition of the gluon fragmentation function 


at high energy colliders from Erst principles. 

V. CONCLUSIONS 

Since the dehnition of the gluon to hadrons fragmentation amplitude < H,X\g > involves 
gluon held Q^°‘{x) corresponding to single gluon incoming state \g >, the present dehnition 
of the gluon fragmentation function at high energy colliders [which involves non-abelian held 
tensor F^'"°‘{x) = — d''Q^°'{x) + gf°‘^^Q^’^{x)Q''^{x) instead of Q^°‘{x) in the initial 

state] is not consistent with the single gluon incoming state \g >■ In this paper we have 
derived the correct dehnition of the gluon fragmentation function at high energy colliders 
from hrst principles which is obtained from the single gluon incoming state \g > and is gauge 
invariant and is consistent with the factorization theorem in QCD. 
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